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Abstract 

In the statistical analysis of shape a goal beyond the analysis of static 
shapes lies in the quantification of 'same' deformation of different shapes. 
Typically, shape spaces are modelled as Riemannian manifolds on which 
parallel transport along geodesies naturally qualifies as a measure for the 
'similarity' of deformation. Since these spaces are usually defined as com- 
binations of Riemannian immersions and submersions, only for few well 
featured spaces such as spheres or complex projective spaces (which are 
Kendall's spaces for 2D shapes), parallel transport along geodesies can 
be computed explicitly. In this contribution a general numerical method 
to compute parallel transport along geodesies when no explicit formula 
is available is provided. This method is applied to the shape spaces of 
closed 2D contours based on angular direction and to Kendall's spaces of 
shapes of arbitrary dimension. In application to the temporal evolution 
of leaf shape over a growing period, one leaf's shape-growth dynamics can 
be applied to another leaf. For a specific poplar tree investigated it is 
found that leaves of initially and terminally different shape evolve rather 
parallel, i.e. with comparable dynamics. 
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1 Introduction 

For more than two millenia, the analysis of form of biological entities has been an 
enticing object of human occupation. While early work tended to be speculative 
in nature, the current state of mathematics and computational power allow 
to develop and simultaneously verify theoretical results, thereby increasingly 
driving scientific progress as witnessed today. 

This present work has been motivated by joint research with the Institute 
for Forest Biometry and Informatics at the University of Gottingen, to compare 
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leaf growth dynamics within single specimen, species and taxa for identification 
of gene expression. The endeavor is challenging as it touches problems at least 
as old as The o phras tus' (371 - 287 b.C.) famous book on "plant growth", cf. 
Theophrastul (jl976l) . 



The first step of this project is the subject of this work: to develop a frame- 
work allowing to compare shape dynamics. To this end, we model biologi- 
cal growth by (generalized) geodesies in shape space. We do so because the 
geometry of shape spaces in which travel along geodesies requires no energy 
seems linked to the physiological reality of growth pref erring to minimize en - 



ergy. This "geodesic hypothesis" , originally stated by iLe an d Kume 



is further suppor ted by earlier research, cf. iHuckemann and Ziezold 2006); 



(2000)) 



Hotz et all ( 2010l) . As with the "geodesic hypothesis" one can say that geodesic 



shape deformation of two different shapes is the "same" if the impetus of the 
first deformation is transplanted to the second with no loss of energy. In the lan- 
guage of Riemannian geometry this translates to the condition that the initial 
velocity of the second geodesies is the parallel transport of the initial velocity of 
the first geodesic. If the deformations are not the same, i.e. the geodesies are not 
parallel at the first and the second shape, this concept gives a correlation-based 
distance between the deformations. 

In consequence, the aim of this paper is to provide for parallel transport on 
shape spacs. Recall that most shape spaces can be viewed as Riemannian im- 
mersions or submersions or, combinations thereof. Explicit formulae for parallel 
transport are only availabl e for special space s , e.g. for spheres and Kendall's 
spaces of planar shapes, cf. Huckemann et al. ( 20091 ). In general, parallel trans- 



port may be difficult to compute and be only available numerically. In the 
following Section [2] we provide for a general method to compute parallel trans- 
port on shape spaces. In view of our application the method is illustrated 
in Section [3] for the spaces of closed 2D co ntours based on angular direction 
wi th and without spec ific initial point (cf. IZahn and Roskiesl (1972) as well 



as 



Klassen et al.l (l2004f)1. and in Sectio n 2] for Kendall's landmark based shape 



spaces (e.g. Drvden and Mardia ( 19981 )) 



In Section [5j we compare parallel transport on the spaces of closed contours 
with parallel transport on Kendall shape spaces for simple regular polygonal con- 
figurations. While all sectional curvatures in Kendall's shape space are bounded 
from below by 1, it turns out that the corresponding subspace of closed contours 
is flat. 

Finally in Section leaf growth of one leaf is transported parallelly to other 
leaves and both shape evolutions are compared with one another. For a spe- 
cific Canadian black poplar investigated we find that leaves with initially and 
terminally different shapes tend to evolve parallel, in particular so if no shape 
anomalies are present. Thus the geodesic hypothesis can be extended to the 
parallel hypothesis: 

biological growth of related objects, possibly of initially and terminally different 
shape, tends to follow parallel geodesies, 

Using Euclidean approximations in landmark based shape spaces rather than 
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geodesies, this hypothesis was originally coined bv lMorris et al. I d2000h who ob- 
served parallel growth patterns. Readers primarily interested in the application 
can directly skip to Section O 



2 Parallel Transport (PT) 

This section begins with a review of basic concepts of Riemannian geometry 
found in any standard textbook (specifically iLand ( 19991 ) is very appropriate for 



the following), in particular formulae relating covariant derivatives of Rieman- 
nian immersions and submersions. These provide differential equations lifting 
the parallel transport on shape space to Euclidean or Hilbert space. 

For a Riemannian manifold M, possibly of countable dimension denote by 
(V p , W p ) M the Riemannian metric of tangent spaces and by Vy W the covariant 
derivative of vector-fields. Here V, W € T(M) denote vector-fields with values 
V p ,W p in the tangent space T p M of M at p € M. d M (p,p') is the induced 
metrical distance on M for p, p' E M, V <E> W denotes the outer product defined 
by (V ®W) X = (X, V) W. A vector-field W G T(M) is parallel along a smooth 
curve t — > j(t) on M if it satisfies the ordinary differential equation (ODE) 

V^W = 0. (1) 

It is well known that there is locally a unique solution W along 7 for a given 
initial value. In Euclidean or Hilbert space the left hand side has the simple 
form ©. 

In particular, geodesies are characterized by the fact that their velocity is 
parallel: V77 = . 

The covariant derivative is often called a covariant connection. Indeed, if 
two offsets p,p' G M can be joined by a unique geodesic segment of minimal 
length, their respective tangent spaces are connected via parallel transport (PT). 

Definition 2.1. w' 6 T P 'M is the parallel transplant of w € T p M if there are 

1. a unique unit speed geodesic t — ¥ j(t) connecting p = 7(0) with p' = 
1 {d M (p 1 p')), and 

2. a vector field W € T(M) parallel along 7 with W p — w, W p > = w' . 

A sufficient condition for the existence of such a unique connecting geodesic 
is that M is finite dimensional and p' is sufficiently close to p. In case of 
infinite dimension, examples of complete spaces can be constructed which do not 
feature minimizing geodesies between arbitrary close points (e.g. Lang ( 19991 



pp. 226/7)). For our applications in mind this fact seems less troublesome since 
infinite dimensional spaces considered here are built from projective limits of 
finite dimensional spaces. 

The Euclidean and Hilbert spaces K™ (for Hilbert space n = 00) can be 
identified with all of their tangent spaces, i.e. (v,w) R = X)"=i an( ^ ^ ne 
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covariant derivative is just the usual multivariate derivative by components, 

vr; 1 ,..^ ) (- i ,---,-")=E^(^---^J- 

In particular, if v — i.e. = v 1 we have that 

nl)W = jW x(t) . (2) 

Thus as desired, parallel transport on Euclidean and Hilbert spaces is given 
by afhne translations. 

For short we write W(t) for the value of W along a selfunderstood smooth 
curve t — > 7(f) and 

W(t) := d -W l(t) 

in the Euclidean/Hilbert case. 

A surjective linear mapping / : E — >• F of topological vector-spaces spZifs 
in F if kern(/) has a closed complement F in E such that F x kcrn(/) = E 
as topological vector-spaces, in particular, kern(J) — »■ F — > F is a s/iort ea;ac< 
sequence. Another wording is that F splits over E. 

A smooth mapping $ : M — > A of Riemannian manifolds M and A induces 
a differential mapping d$ p : T p M — > T^^N of tangent spaces. $ is called 

an immersion if i> is injective and if every T$( p ) A splits over d& p T p M, 

a submersion if $ is surjective and if every d$ p splits in T$( p )A, 

an isometry if (Vp,Wp) M = (d$ p V p , d^lf/ , Vp G M and V, W e 
T(M) . 

An isometric immersion (submersion) is a Riemannian immersion (submersion) 
respectively. 

Riemannian Immersions. If $ : M — > N is a Riemannian immersion then 
the tangent spaces of A split into the tangent spaces of $(M) and its orthogonal 
complements, the normal spaces 

I*(p)JV = F $(p) $(M)® A„ (p) $(M). 

As a consequence of the implicit function theorem, every Riemannian immersion 
$ : M -> A admits locally an implicit representation $:[/flJV4 A $ ( p )$(M) 
such that nM) = [/fl $(A). Here [/ is a suitable neighborhood of 
in A. Hence, we have with A, Y e T(M) and arbitrary local extensions A, Y e 
T(t/ n A) of d$A, d$Y e TM that 

(idr(iv) -d*) (Vf Y) = d$(V^Y). (3) 
In particular, d^$( p ) spans the normal space A$( p )<I>(M). 
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Theorem 2.2. Suppose that an embedding idu ■ M M" is a Riemannian 
immersion in Euclidean (n < oo) or Hilbert space (n = oo), t — > ^(t) a geodesic 
in M , {Vj(t) : j G J} an orthonormal smooth base for N 1 ^M and W a vector- 
field in M. Then W is parallel along 7 if and only if it satisfies the linear 
differential eguation 

W(t) = - W[t). 

Proof. The assertion is an immediate consequence of ([3]) and the fact that 

= ^(W(t),Vj{t)) = {W{t),V ] {t)) + (W{t),V ] {t)) 
for all j € J by hypothesis. □ 

Riemannian Submersions For a Riemannian submersion $ : M — > Q from 
the top space M to the bottom space Q, tangent spaces split as follows: every 
fiber <I> _1 ((7), q G Q is a submanifold of M that is locally a topological embed- 
ding. With the vertical space T p <£> -1 (<£>(p)) along the fiber and its orthogonal 
complement, the horizontal space, we have 

T p M = Tp®- 1 ($(p)) © H P M . 

Since H p M = T$r p \Q, every V G T(Q) has a unique horizontal lift V G H P M 
characterized by d$y = V. For arbitrary PF G T(M) denote by : p — > Wp 1 
the orthogonal projection to the v ertical space. 

The following Theorem due to lO'NeiUl (|l966tl (cf. also|Lang| (|l999l p.386)) 



allows to lift bottom space parallel transport to the top space. In addition to 
([3]) this provides the vertical (normal) part as well, which is in general non-zero 
for submersions. 

Theorem 2.3. Let <I> : M — > Q be a Riemannian submersion and let X,Y G 
T{Q). Then we have with the Lie bracket [-, •] on M that 

vfy = v|y+ \ [X,Y] X . 

We are now ready for the ODE of parallel transport on a Riemannian im- 
mersion followed by a Riemannian submersion. 

Theorem 2.4. Suppose that $1 : M R™ is a Riemannian immersion in 
Euclidean (n < 00 ) or Hilbert space (n = 00), : M —> Q a Riemannian 
submersion and let W be a vector field on M horizontal along a horizontal 
geodesic j(t) on M . Then d&2W is parallel along $2 j{t) if and only if 

w(t) = - w(t)-j2ck>i(j(t),w(t)). 
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Here, {Vj(t) : j € J} denote an orthonormal smooth base for the normal space 
iV 7 ( t )M C K™ and duil are the exterior derivatives of an orthonormal and smooth 
base {Uk(t) : k € A^} of the vertical space Ty^[j(t)} C T 7 ( t )M for suitable index 
sets J and K . 

Proof. Suppose that we have a vector field X 6 T(Q) with horizontal lift X(t) = 
■y(t). If $2 ° 7 is geodesic and c^VF parallel with horizontal lift W, Theorem 
O yields 

= ^Z)<[^.W].^*>]^* = -5Z^fc(7,W) (4) 

keK k£K 

making use of the well known (e.g. Lane ( 19991 . p. 126/7)) 



([X,W],U k ) =X(U k ,W)-W(U k ,X) -2doj k (X,W) 

with the exterior derivative duk of the one- form Wk dual to Uk- On the other 
hand, since d$i : T P M -> T P M C M" is the identity, formula © yields 

V^W M = W{t)- i^Vj{t)®V {t^ W{t) 

as in the proof of Theorem l2.2l . Putting together (|4]) and ([5]) gives the assertion 
of the Theorem. □ 



3 PT for Closed 2D Contours 



We define the two shape spaces of closed 2D constant - speed contours based on 
angular direct ion as introdu c ed by Zahn and RoskiesI (1972) in the geometric 
formulation of Klassen et al.l ( 2004 ) . 

Suppose that z : [0, 2tt] — > C, s t— > z(s) is a constant-speed parameterization 
of a smooth, closed, curve of length L winding once counterclockwise around 
each interior point. Let 



6(s) 



rg (/(«)) -arg(z'(0)) 



with 



J( s ) = J± e »(e(s)+arg(2:'(0))+ S ) 
27T 



(6) 



Obviously, the Zahn-Roskies shape (ZR-shape) 6 is invariant under translation, 
scaling and rotation z(s) — > c + \e^z(s). Moreover, subtracting s (the curves 
to be modelled wind once around their interior) norms 9 such that it is 2n- 
pcriodic. Vice versa, from every converging Fourier series an a.e. differentiable 
constant-speed 2D curve can be reconstructed by integrating ([6]). This curve is 



6 



unique modulo translation, scaling and rotation. Thus a linear subspace of the 
Hilbert space £ of Fourier series is the ZR-pre-shape space 

oo 

Szr ■= \0(s) = ^2 ( Xn cos ( ns ) + 2Msin(ns)) : 

1 X 

n=l 

00 

^x n , y = 0} . 



T 2 



Xo 



71 = 1 



As usual, 2tt(8, 77) := f Q v 9(s)rj(s) ds and ||0|j 2 := (6>, 0). The tangent spaces 

TgSzR are identified with S'z.r C £ 2 . Since the curves in question are closed, 
we have with the non-linear mapping 

«2 „ 



* : 



1 ^ 



Jo e 



(0( S )+s) 



that the ZR-shape space is the implicit sub-manifold 
Szr := € Sztf : #(0) - 0} . 

Obviously, the ZR-shapes of closed not self-intersecting contours form an open 
subset containing the origin, which corresponds to the shape of the circle. 

Additionally considering closed curves invariant under change of initial point 
z{s) — > z(s + so) (e.g. amorphous curves with no preassigned initial point) by 
defining this action of the unit circle S 1 3 sq on Yizr the invariant ZR-shape 
space 



n ZR 



(s^Uo})/^ 



(7) 



is obtained. 

Since rotation and parameter shift are equivalent for circles, the correspond- 
ing invariant ZR-shape is thus a singularity of in fact its only singularity. 

Parallel Transport on Szr and TJ ZR Geodesies on Szr as well as on S J ZR 
betwee n two given p o ints c an be comp uted via a technique c alled geodesic shoot- 
ing ', cf. I Miller et al.l ([20061 ) as well as iKlassen et al.l ([20041 ) , or much faster via 
a variational approach ISchmidt et al. (|2006l) . Since Szr £ is a Rieman- 
nian immersion with the global implicit definition = we have that the 
normal space in Szr at 9 € 'Szr is spanned by Vi(9) = «4 cos (9(s) + s) and 
^2(0) = s h-> sin (6(s) + s). Orthogonalization yields the base 



Wi := 



__L 

|Vil 



W 2 := 



As a consequence of Theorem 



V 2 
we have 



V2-(V2,W 1 )W 1 
{V 2 ,Wi)Wi\ 
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Theorem 3.1. A vector-field WCt) in Tizr is parallel along a geodesic 7 in 
£ zr if and only if it satisfies the linear differential equation 

W(t) = (8) 

In practice, © can be solved numerically by orthogonally projecting to 
T 7 ( t )E2K in every iteration step. 

We now turn to the submersion 0. The vertical space at 

00 

6(s) = x Q + 2J (x n cos(ns) + y n sin(ns)) G Ezr 
n=l 

is spanned (if convergent) by the single vertical unit length direction 



?'(s) D^Li n { - x nd Vn + y n d x 



The exterior derivative of its dual is hence 



V2. 



IDti^ti' yn x n' nn ( n dx n f\dx n — Tidy 71 Ady n ) 

Yin n' nn ( n VnVn' +"-'^ti^ ti ' )dy n /\dx n 



(9) 



In conjunction with Theorem 12.21 Theorem 12.41 and Theorem 13.11 one obtains 
after a tedious computation 

Theorem 3.2. The vector-field 

00 

W(t) = u Q (t)d X0 + ^ K(«)9, B + v n (t)d Vn ) 

71=1 

is a horizontal lift to the top space Ezr. of the bottom space parallel transport 
along a geodesic in E I ZR 

00 

7 S (£) = x (t) + ^ (x n (t) cos(ns) + y n (t) sm(ns)) 

n=l 

horizontal in Ezr if and only if it satisfies the linear differential equation 
W(t) = 

-{j t W ^°^) ®Wi( 7s (t)) +j t W 2 (is(t)) ®W 2 (7.(*))) W(*) 
(<Tl(*)®'>4'(*).W(t)®7.(t)> 

-<Ti(*)»7.(*),W r (*)8'yi'(*)>) 



2||7iW 



r <7i(t),w(t))- 



IWWII 2 
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with the derivatives defined as 

7«(*) = XT=i n ( ~ x n{t)d Vn + y n (t)d Xn ) , 

j s (t) = x {t)d Xo + Y,n=i ( i n(t)d Xn + y n {t)d Vn ) , 

■YM = Y,n=l n i - i n(t)d Vn + i)n{t)d Xn ) , 

<{t) = -En=in 2 (x n (t)d Xn +y n (t)d yn ), 

if convergent, and the inner product defined by 

(Ei ® Ej,E k ® Ei) :— S(i,j),(k,i) 

for an orthogonal system Ej and index set J 3 j. 

In practice, convergence of the series for the derivates is not an issue since as 
remarked earlier, computations are carried out using only finitely many Fourier 
coefficients. 



4 PT for Kendall's Shape Spaces 

Kendall's landmark based similarity shape analysis is based on configurations 
consisting of k > m + 1 labelled vertices in E m called landmarks that do not all 
coincide. A configuration 

x = (x 1 , . . . ,x k ) = (x %J )i<i<m,i<j<k 

is thus an element of the space M(m, k) of matrices with k columns, each an Tri- 
dimensional landmark vector. Disregarding center and size, these configurations 
are mapped to the pre-shape sphere 

M = S k m := {p G M(m, k-1): \\p\\ = 1} , 

where ||p|| 2 = (p,p) and (p,v) := tr(pv T ) is the standard Eucli dean product. 



This c an be done by, say, multiplying by a sub-Helmert matrix, cf. iDrvden and Mardia 
(1998) for a detailed discussion of this and other normalization methods. The 
canonical Riemannian immersion c — > M(m, k — 1) comes with a global im- 
plicit definition ty(x) — \\x\\ — 1 = giving rise to a single normal field. Hence, 
using ([2]) and ((3|), the covariant derivatives along a curve 7 on relate as 

v^ t) w(t) = w(t)-(w(t)Mt))^t) (10) 

As a consequence of (fTO]) and (JXJ) , unit-speed geodesies on are great circles 
of form 7(i) = x cos t + v sin t with x, v £ and (x, v) = 0. 

In order to filter out rotation information define the regular part (S^)* := 
{x £ S k : rank(x) > m — 2} (an open dense subset of S^J and a smooth and free 
action of SO(m) by the usual matrix multiplication (S*J* A (S k n )* : p ^ gp 
for g £ SO(m). The orbit 7r(p) = {gp: g £ SO(m)} is the Kendall shape of 
p £ S^j and the quotient 

tt : (S k n )* -+ (E* )* := {S k m f / SO{m) (11) 
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is called Kendall's shape space. Note that projecting from the entire pre-shape 
sphere S!^ would have led to a non-manifold quotient (m > 3), which is usually 
called Kendall's shape space. We consider here only the regular part such that 
is a Riemannian submersion. With the orthogonal decomposition g[(m) = 
o(m) © SM(m) of the Lie algebra gt(m) = M(m,m), the Lie algebra o(m) of 
skew-symmetric matrices in gl(m) and the vector-space of symmetric matrices 
SM(m) in flt( m) we have the f ollowi ng orthogonal tangent space decomposition 
for ieSj„ cf. iKendall et all (Il999l p. 109). 



gl(m) = o(m) © SM(m) 

_ (12) 



i -x t -x T 



T X S^®N X S^ = T x n(x) © H x S k m ®N x S k m 

For x € (S^)* both mappings are surjective and H x S k n = T 7r ( 2 ,)(I]J;J*. In order 
to compute the horizontal lift of bottom space parallel transport as in Theorem 
3.21 we need an orthonormal base for the (m(m — l)/2)-dimensional vertical 
space T x ir(x) and the exterior derivative of its duals. From (|12|) we have at once 
a (in general not-orthogonal) base {eijX : 1 < i < j < m} with base system 



1 for a = i, j3 = j 
with e a ^ = { -1 for a=j,/3 = i 
else 



of o(m). From the former obtain an o.g. base system {Vij(x) : 1 < i < j < m} of 
T x tt(x) through Gram-Schmidt orthogonalization and let tOij(x) be the one-form 
dual to Vij(x) (1 < i < j < m). For m = 2 there is a single vertical unit-direction 
V"i2(cc). For m > 2, however, duji2(x) is about as complicated as ([9]), and even 
for m = 3, the other two derivatives and their application to vector-fields result 
in expressions to o leng thy to be written down, cf. also the rather complicated 



examples in |LgJ (|2003l ). Using a computer algebra program, however, these 



expressions and their respective values can be easily computed by symbolic 
differentiation. Hence, Theorem 12.41 and (fTOl) yield th e follow ing Theorem. The 



special case m = 2 is taken from Huckemann et all ( 20091 Theorem A. 6), cf. 



also|LgJ (|2003l Theorem 2) 



Theorem 4.1. A vector-field W{t) is a horizontal lift to the top space of the 
bottom space (Sj^)* parallel transport along a geodesic ^{t) — a; cos £ + v sin t 
horizontal in (S^)* if and only if it satisfies the ODE 

w(t) = (w(t)Mt)) 7(*) 

- d ^ (?(*)) (">(*)' W(tj) Vij ( 7 (*)) . 

l<i<j<m 

For m — 2 this ODE has the explicit solution 

W(t) = W(0)- ((W(Q),v)v + (W(0),e 12 v)ei2V 
+(W(0), v) % v (t) + (W(0), e 12 v) ei,37*,»(*) 
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(c) From C2 to the parallel trans- 
plant of o~3 from Figure \l(a)\ 



(c) From 02 to the parallel trans- 
plant of 0-3 from Figure \2(a)\ 



Figure 1: Equidistant deformation 
along geodesies in S zr 2 ■ The bul- 
let marks the pre-assigned initial 
point. Figure 1(c) depicts the par- 



allel transplant to a 2 of" the geodesic 



from Figure 1(a) along the geodesic 
depicted in Figure 1(b) 



Figure 2: Equidistant deformation 
along geodesies in Kendall's land- 
mark based shape space S|. The 
kinks signify landmarks. Notation 
as in Figured 



5 Curves with Rotational Symmetry 

Denote by £ ZR k the sub-space of c losed curves with fc - fold r otational symmetry. 
The following is an observation of Zahn and Roskiesl (119721) : 



Theorem 5.1. 9 represents a closed curve with k-fold rotational symmetry 
k > 1 if and only if x n ,y n = for all n ^ mod k. 
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Since the arithmetic sum of closed curves with /c-fold rotational symmetry 
is again of fc-fold rotational symmetry we have at once: 

Corollary 5.2. For each k = 2,3 . . ., ^zR k is a fl a t linear submanifold of T,zr 
and parallel transport on T> ZR k is affine. 

Within the subspace T, ZR 2 of closed curves with two-fold symmetry consider 
a simple example of three shapes: two rectangles o~i and o~% differing only by 
their initial point and a hexagon 03. The 'same' (i.e. parallel) deformation from 
(Ti to CT3 is applied to 02- 

Figure [1] illustrates that deformation in the geometry of Zahn-Roskies' shape 
space, Figure [2] gives it in the geometry of Kendall's shape space. Comparing the 
respective Subfigures (a) and (b) over the different geometries shows that the 
geodesic deformation with fixed initial and terminal shape gives almost identical 
intermediate shapes. 

In the respective Subfigures (c), the difference between 'same' shape de- 
formation over the two geometries is hardly notable in the beginning of the 
deformation. Near the end, however, it becomes notable: number the kinks 
(the landmarks in E^) counterclockwise from 1 (bullet) to 6 and denote by 
the line connecting the i-th kink with the j-th kink. Then (2, 5) remains par- 
allel to (6, 1) and (3,4) in the Kendall geometry whereas in the Zahn-Roskies 
geometry it turns in direction beyond (2, 3) and (5, 6). 



6 Parallel Leaf Growth 



Let us quickly overview some very recent developments within two millenia of 
research on plant form. With the application below in mind, we are interested 
in a flexible and realistic representation of leaf contour shape. Flexibility in this 
context means that we are looking for a model in which nature not only chooses 
values of parameters in a pre-defined parameter space but rather the parameter 
space itself. We caution that such a non-parametric model may come at the 
cost that nature's parameters may not be simply geometrically interpretable 

Parametric in this sense are the we ll estab li shed models in volving allometry, 
still of interes t toda y: e.g. Gurevitch ( 19921 ): iBurtonl ( 20041 ). or the superfor- 
mula of Giclis (2003). Al so, models based on landmarks such as lDickinson et all 



(1987); 



Jensen! (1990) or I Jensen et al.1 ( 20021 ). can be viewed as projecting na- 



ture to a pre-specified parameter space by leaving out the parts of the contour 
between landmarks. Note that parametric models are highly successful e.g. for 
lant classification, genetic hybrid identification, cf. IJorgensen and Mauricio 
2005), or in the Climate Leaf Analysis Multivariate Program (CLAMP) of 
Wolfd ([19931 ) whi ch is fundamental to paleoclimate and present day climate 
reconstruction, cf. lEndress et al. I (120001) . 

On the other hand, models building on the shape spaces of Zahn and Roskies 
(cf. Section 13]) are non-parametric, even though they are not entirely free of con- 
straints: in view of landmark-based shape analysis (cf. Section 0]), restricting 
to unit speed velocities translates into infinitesimally placed landmarks. The 



12 



different geometry, however, liberates from the necessity to identify homologous 
landmarks, by imposing infinitesimal uniform growth. The latter is certainly 
debatable. Curiously, such non - param etric models introduced as eigenshape 
analysis (building on iLohmaiml ( 19831) ) have initially stirred controversy be- 
cause parameters were not simply geometrically interpreta ble an d beca use with 
lacking initial point, registration was not satisfactory, cf. Rohlj f|l986h . While 
the former is precisesly a desired feature , introducing the geo metric concept of 
the quotient YJ ZR = (E Z r \ {O})/^ 1 bv lKlassen et al.l f|2004h settles the latter 
objection. It seems, however, that the natur al non-Eucl i dean geometry of Hzr 
is not fully r ealized in the community, cf. 
iHearnl (120091) . 



Ravi (|1992I ); iKrieger et al.l (|2007l ) 



For sake of completeness, even though not practicable for our purpose be- 
cause of high sensitivity to boundary noise, let us briefly mention a third ap- 
proach of shape modeling based on the leaf' s vein structure. W ith methods for 
automated venation extraction available (cf. iFu and Chi ( 20061 ) ), although com - 



(2009) 



1978) 



putationally much more challenging than contour ext raction. ILu et al 
link vein structure to the concept of shape spaces b y iBlum and Nagel 
based on medial skeletons, cf. also Pizer et al. ( 20031 ). Undoubtedly, modeling 
the vein structure gives deep insight into physiological, hydraulical and biomc- 
chanical aspects of leaf formation. Parameter spaces thus obtained should be 
closest to nature in the above sense. Current research in venation patterns, 
however, shows that l eaf shape diversification is still poorly understood (e.g. 



Niinemets et al.l (|2007)) 



Obviously, for our purpose of modeling entire leaf contours while being as 
non-parametric as possible, the space Tizr suits ideally. For the problem at 
hand there is no need for pre-registration as the leaves in question are naturally 
aligned by petiole (the base point where the stalk enters the blade forming the 
main leaf vein) and apex (the terminal point of the main vein, usually the leaf 
tip) location. One could almost equivalently align by petiole location and the 
initial direction of the main leaf vein. 



be1b8: original 




be1b8: projection to its geodesic 




Figure 3: Shape evolution of a black poplar leaf over two weeks. Left: original 
contours. Right: contours obtained from projecting to geodesic evolution. 

In an application we consider five leaves of a Canadian black poplar tree at an 
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Figure 4: Shape evolution of black poplar leaves over two weeks. Left column: 
original contours. Right column: contours obtained from traversing the first 
contour along the parallel translate of the geodesic from Figure [3J 
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experimental site of the Department of Forest Biometry and Tree Physiology 
of the University of Gottingen. Their contours have been non-destructively 
extracted over their growing period at 10 approximately evenly spaced days. 
For the following computations we have considered the Euclidean (2N + 1)- 
dimensional subspace of £ 2 using Fourier coefficients, (xq, x±, yi, . . . , xn, Vn) 
with N = 100. 

Leaf "belb8" (left image of Figure[3]) exhibits the most regular sh ape. In con- 
cord w ith earlier observations of different leaves using landmarks (cf. iHotz et all 
(2010)). the temporal evolution occurs almost along the geodesic determined by 
initial and terminal shape (as depicted in the right image of Figure [3]). The 
initial direction of this geodesic has been parallely transplanted to the initial 
shapes of leaves "belb9" , "belb3" , "be2b4" and "be2bl" . In the right columns 
of Figure 2J the shapes along these new geodesies starting at the corresponding 
initial shapes have been recorded at the corresponding points in time. The left 
columns of these figures depict the original temporal shape evolution. 

The common shape dynamics displayed by the original leaf contours (left 
columns of Figures [3] and [4} seems two-fold. First, an increase of base angle. 
Second, different growth ratios are not visible at the apex as its angle remains 
nearly unchanged. Individual effects are non-symmetric and non uniform lateral 
growth. Also, leaf 'belb3' develops a notch left, slightly below the apex, for leaf 
'bc2bl' an original notch also left, slightly below the apex attenuates. 

Obviously (right columns of Figures [3] and [4]) leaves "belb9" , "belb3' and 
l be2b4" follow rather closely the parallel transplant of the geodesic of leaf 
"belb8" . Original non-uniform growth is uniformized and, stronger than origi- 
nally, apexes acuminate. Even though all of their initial and terminal shapes are 
quite different, one can say that their temporal evolution is rather similar. This 
seems to be less the case for leaf "be2bl" . Its observed growth tends to eliminate 
its initial strong dent at north-west-north while along the transplanted geodesic, 
this dent remains, causing increased distal growth at the tip. One could argue 
that in order to restore an original contour defect, natural growth deviates from 
its "original" plan. Certainly, such phenomena deserve future research. 



b 


bclb8 belb3 be2b4 be2bl 


P(v h , Wb) 
fJ-(vb, Wb) 


0.17 0.12 0.44 0.083 
0.99 0.96 1.0 0.88 



Table 1: Measuring par allelity of geodesies at first shape. Top row: non-central 
correlation |1^[ ) of Fourier coefficients of Vb ( initial velocity of geodesic approx- 
imating leaf shapes of leaf 'b' at its first shape) and parallel transplant Wb of 
fbebis (initial velocity of geodesic approximating leaf shapes of leaf 'belb8' at its 
first shape) to the first shape of leaf 'b'. Second row: the same for the (1 — p)- 
values obtained from \13\) . 

As a measure for parallelity, the cosine of the angle between initial velocity of 
geodesies or equivalently the the correlation of the respective Fourier coefficients 
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could be taken, cf. first row in Table [T] At first glance, in contrast to Figure 21 
these numbers suggest an almost non-existent amount of parallelity. On closer 
inspection, taking into account, however, that the vectors compared are high- 
dimensional (of dimension 2N + 1 = 201), the first and third number of the 
second row in Table [1] indicate high correlation as expressed in the third row: 
if two random vectors v, w would be independently sampled from a uniform 
distribution on the (n — l)-dimensional unit-sphere, then the density of their 
angle is proportional to the surface of the (n — 2)-hypersphere with radius sin <fi 
determined by this angle. Based on this consideration we propose the following- 
measure for parallelity of random v, w G R n 

J v ^ ' sin" z <j> d(j) 
H(v,w) = 1 n , n _ 2 , with (13) 



r 7T . fi 

Jo sin 

P{v,w) = ■ 14) 

\\v\\ \\w\\ 

Finally, let us note that the curves of the contours of the leaves be2b4 and 
be2bl projected along the transplanted geodesic of leaf belb8 start to self- 
intersect as the leaves grow older. While the effect is rather small, it may 
indicate that the natural metric of Y*zr be adjusted in order to maintain the 
hypothesis of geodesic growth. 



7 Discussion and Outlook 

In this exposition a method to compare shape dynamics has been proposed based 
on parallel transport of geodesies. While there is quite a few work available, 
modell ing temporal shape evo l ution by spe c ific cu r ves and spline s in sh ape space 
((e.g. IJupp and Kentl (|l987l ); iKent et all (|200ll) : iKume et alj (|2007h ). to the 



knowledge of the author this is the first time that dynamical aspects of different 
shapes have undergone a comparison based on the intrinsic geometry of shape 
space. 

In application to Botany, growth of leaves of different shapes has been com- 
pared. As underlying shape representation the space Y^zr of closed contours 
based on angular direction has been employed. Within this space, in contrast 
to other models, entire leaf contours can be retrieved in a non-parametric way. 
In a simple toy example, parallel transport on Yizr seemed locally similar to 
parallel transport on Kendall's shape spaces. 

This similarity can be rephrased as saying that the landmarks have been 
"correctly chosen". Parallel transport thus may serve as a tool to address an 
open problem in landmark based shape analysis: optimize number and location 
of landmark placement for a specific problem at hand. Obviously, too few and 
wrongly placed landmarks have low predictive power, while too many landmarks 
reduce power due to undesired variation. Dealing with this latter effect usually 
requires further methodology, e.g. statistical regularization. 
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In the context of dynamical shape analysis for leaves, under optimal land- 
mark placement, since computation of geodesies, e.g. geodesic PCA, is compu- 
tationally much faster on Kendall's shape spaces than on Sz_r, one can perform 
parallel transport on Kendall's shape spaces and obtain the complete bounding 
contour in a non-parametric way by mapping to Tizr- This lays out a path for 
the future steps of the challenging endeveavor of statistically comparing shape 
dynamics laid out in the Introduction. 
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